Scalar const c0 = pow(a, 2);
Scalar const c1 = pow(b, 2);
Scalar const c2 = c0 + c1;
Scalar const c3 = pow(c2, 5.0/2.0);
Scalar const c4 = c2 + 1;
Scalar const c5 = atan(sqrt(c2));
Scalar const c6 = pow(c5, 2);
Scalar const c7 = c6*k[0];
Scalar const c8 = pow(c5, 4);
Scalar const c9 = c8*k[1];
Scalar const c10 = pow(c5, 6);
Scalar const c11 = c10*k[2];
Scalar const c12 = pow(c5, 8)*k[3];
Scalar const c13 = 1.0*c4*c5*(c11 + c12 + c7 + c9 + 1.0);
Scalar const c14 = c13*c3;
Scalar const c15 = pow(c2, 3.0/2.0);
Scalar const c16 = c13*c15;
Scalar const c17 = 1.0*c11 + 1.0*c12 + 2.0*c6*(4*c10*k[3] + 2*c6*k[1] + 3*c8*k[2] + k[0]) + 1.0*c7 + 1.0*c9 + 1.0;
Scalar const c18 = c17*pow(c2, 2);
Scalar const c19 = 1.0/c4;
Scalar const c20 = c19/pow(c2, 3);
Scalar const c21 = a*b*c19*(-c13*c2 + c15*c17)/c3;
result[0] = c20*f[0]*(-c0*c16 + c0*c18 + c14);
result[1] = c21*f[0];
result[2] = c21*f[1];
result[3] = c20*f[1]*(-c1*c16 + c1*c18 + c14);
